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Abstract 

We provide a number of new construction techniques for cubical complexes and cubical 
polytopes, and thus for cubifications (hexahedral mesh generation). As an application we 
obtain an instance of a cubical 4-polytope that has a non-orientable dual manifold (a Klein 
bottle). This confirms an existence conjecture of Hetyei (1995). 

More systematically, we prove that every normal crossing codimension one immersion of 
a compact 2-manifold into is PL-equivalent to a dual manifold immersion of a cubical 
4-polytope. As an instance we obtain a cubical 4-polytope with a cubation of Boy's surface 
as a dual manifold immersion, and with an odd number of facets. Our explicit example has 
17 718 vertices and 16 533 facets. Thus we get a parity changing operation for 3-dimensional 
cubical complexes (hexa meshes); this solves problems of Eppstein, Thurston, and others. 
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1 Introduction 

A d-polytope is cubical if all its proper faces are combinatorial cubes, that is, if each k-face of 
the polytope, k € {0, . . . , d — 1} is combinatorially equivalent to the /c-dimensional standard 
cube. 

It has been observed by Stanley, MacPherson, and others (cf. [3] [20]) that every cubical d- 
polytope P determines a PL immersion of an abstract cubical (d — 2)-manifold into the polytope 
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boundary dP = S'^~^. The immersed manifold is orientable if and only if the 2-skeleton of the 
cubical d-polytope {d > 3) is "edge orientable" in the sense of Hetyei, who conjectured that there 
are cubical 4-polytopes that are not edge-orientable [17, Conj. 2]. 

In the more general setting of cubical PL (d — l)-spheres, Babson and Chan [3] have observed 
that every type of normal crossing PL immersion of a (d — 2)-manifold into a {d — l)-sphere 
appears among the dual manifolds of some cubical PL (d — l)-sphere. 

No similarly general result is available for cubical polytopes. The reason for this may be 
traced/blamed to a lack of flexible construction techniques for cubical polytopes, and more 
generally, for cubical complexes (such as the "hexahedral meshes" that are of great interest in 
CAD and in Numerical Analysis). 

In this paper, we develop a number of new and improved construction techniques for cubical 
polytopes. We try to demonstrate that it always pays off to carry along convex lifting functions of 
high symmetry. The most complicated and subtle element of our constructions is the "generalized 
regular Hexhoop" of Section 6.4, which yields a cubification of a d-polytope with a hyperplane 
of symmetry, where a (suitable) lifting function may be specified on the boundary. 

Our work is extended by the first author in [28], where additional construction techniques for 
cubifications (i. e. cubical subdivisions of d-polytopes with prescribed boundary subdivisions) 
are discussed. 

Using the constructions developed here, we achieve the following constructions and results: 

• A rather simple construction yields a cubical 4-polytope (with 72 vertices and 62 facets) for 
which the immersed dual 2-manifold is not orientable: One of its components is a Klein bottle. 
Apparently this is the first example of a cubical polytope with a non-orientable dual manifold. 
Its existence confirms a conjecture of Hetyei (Section 5). 

• More generally, all PL-types of normal crossing immersions of 2-manifolds appear as dual 
manifolds in the boundary complexes of cubical 4-polytopes (Section 7). In the case of non- 
orientable 2-manifolds of odd genus, this yields cubical 4-polytopes with an odd number of 
facets. From this, we also obtain a complete characterization of the lattice of /-vectors of 
cubical 4-polytopes (Section 9). 

• In particular, we construct an explicit example with 17 718 vertices and 16 533 facets of a 
cubical 4-polytope which has a cubation of Boy's surface (projective plane with exactly one 
triple point) as a dual manifold immersion (Section 8). 

• Via Schlegel diagrams, this implies that every 3-cube has a cubical subdivision into an even 
number of cubes that does not subdivide the boundary complex. Thus for every cubification 
of a 3-dimensional domain there is also a cubification of the opposite parity (Section 10). This 
answers questions by Bern, Eppstein, Erickson, and Thurston [5] [10] [31]. 

2 Basics 

For the following we assume that the readers are familiar with the basic combinatorics and 
geometry of convex polytopes. In particular, we will be dealing with cubical polytopes (see 
Griinbaum [15, Sect. 4.6]), polytopal (e.g. cubical) complexes, regular subdivisions (see Ziegler 
[33, Sect. 5.1]), and Schlegel diagrams [15, Sect. 3.3] [33, Sect. 5.2]. For cell complexes, barycen- 
tric subdivision and related notions we refer to Munkres [24]. Suitable references for the basic 
concepts about PL manifolds, embeddings and (normal crossing) immersions include Hudson 
[19] and Rourke & Sanderson [27]. 
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2.1 Almost cubical polytopes 



All proper faces of a cubical d-polytope have to be combinatorial cubes. We define an almost 
cubical d-polytope as a pair {P,F), where F is a specified facet of P such that all facets of P 
other than F are required to be combinatorial cubes. Thus, F need not be a cube, but it will 
be cubical. 

By C(P) we denote the polytopal complex given by a polytope P and all its faces. By C{dP) 
we denote the boundary complex of P, consisting of all proper faces of P. If P is a cubical 
polytope, then C{dP) is a cubical complex. If (P, F) is almost cubical, then the Schlegel complex 
C{dP)\{F} is a cubical complex that is combinatorially isomorphic to the Schlegel diagram 
Schlegel(P, F) of P based on F. 



2.2 Cubifications 

A cubification of a cubical PL {d — l)-sphere S'^~^ is a cubical (i-ball B'^ with boundary S'^~^. 
A double counting argument shows that every cubical (d — l)-sphere that admits a cubification 
has an even number of facets. Whether this condition is sufficient is a challenging open problem, 
even for d = 3 (compare [5], [10]). 



2.3 Dual manifolds 



For every (pure) cubical d-dimensional complex C, d > 1, the derivative complex is an abstract 
cubical cell (d — l)-dimensional complex !){€) whose vertices may be identified with the edge 
midpoints of the complex, while the facets "separate the opposite facets of a facet of C," that 
is, they correspond to pairs (F, [e]), where F is a facet of C and [e] denotes a "parallel class" 
of edges of F. This is a cell complex with /i(C) vertices and [d — l)fd-i{C) cubical facets of 
dimension d — 1, d — 1 of them for each facet of C. Hence the derivative complex I'(C) is pure 
{d — l)-dimensional. See Babson & Chan [3, Sect. 4]. 




In the case of cubical PL spheres (for instance boundary complexes of cubical polytopes) , or cu- 
bical PL balls, the derivative complex is a (not necessarily connected) manifold, and we call each 
connected component of the derivative complex 'D(P) of a cubical complex C a dual manifold 
of C If the cubical complex C is a sphere then the dual manifolds of C are manifolds without 
boundary. If C is a ball, then some (possibly all) dual manifolds have non-empty boundary 
components, namely the dual manifolds of dC. 
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The derivative complex, and thus each dual manifold, comes with a canonical immersion into 
the boundary of P. More precisely, the barycentric subdivision of T)[P) has a simplicial map 
to the barycentric subdivision of the boundary complex dP, which is a codimension one normal 
crossing immersion into the simplicial sphere sd(C(9P)). [Normal crossing means that each 
multiple-intersection point is of degree k < d and there is a neighborhood of each multiple 
intersection point that is PL isomorphic to (a neighborhood of) a point which is contained in k 
pairwise perpendicular hyperplanes.) 

Restricted to a dual manifold, this immersion may be an embedding or not. 




Figure 2: The cubical octahedron Og (the only combinatorial type of a cubical 3-polytope with 8 facets), and 
its single immersed dual manifold. 



In the case of cubical 3-polytopes, the derivative complex may consist of one or many 1-spheres. 
For example, for the 3-cube it consists of three 1-spheres, while for the "cubical octahedron" Og 
displayed in Figure 2 the dual manifold is a single immersed (with 8 double points). 

In the case of 4-polytopes, the dual manifolds are surfaces (compact 2-manifolds without bound- 
ary). As an example, we here display a Schlegel diagram of a "neighborly cubical" 4-polytope 
(with the graph of the 5-cube), with /-vector (32, 80, 96, 48). 




Figure 3: A Schlegel diagram of the "neighborly cubical" 4-polytope Cf with the graph of the 5-cube, and its 
dual torus. All other dual manifolds are embedded 2-spheres. 



According to Joswig & Ziegler [21] this may be constructed as 

C| := conv((Q X 2Q) U (2Q X Q)), where Q = [-1, -M]^. 

Here the dual manifolds are four embedded cubical 2-spheres S*^ with /-vector (16, 28, 14) — of 
two different combinatorial types — and one embedded torus T with /-vector (16,32, 16). 
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2.4 Orientability 



Let P be a cubical d-polytope {d > 3). The immersed dual manifolds in its boundary cross the 
edges of the polytope transversally. 

Thus we find that orientability of the dual manifolds is equivalent to the possibility to give 
consistent edge orientations to the edges of the P, that is, in each 2-face of P opposite edges 
should get parallel (rather than antiparallel) orientations; compare Hetyei [17]. Figure 4 shows 
such an edge orientation for a cubical 3-polytope (whose derivative complex consists of three 
circles, so it has 8 consistent edge orientations in total). 




Figure 4: Edge-orientation of (the Schlegel diagram of) a cubical 3-polytope. The edges marked on the right 
must be oriented consistently. 

One can attempt to obtain such edge orientations by moving to from edge to edge across 2-faces. 
The obstruction to this arises if on a path moving from edge to edge across quadrilateral 2-faces 
we return to an already visited edge, with reversed orientation, that is, if we close a cubical 
Mobius strip with parallel inner edges, as displayed in the figure. (Such an immersion is not 
necessarily embedded, that is, some 2-face may be used twice for the Mobius strip.) 

Wo , Wi , W2 Wi_i — , Wi = Vo 




vq — > vi — y V2 vi^i — , ve = wo 

Figure 5: A cubical Mobius strip with parallel inner edges. 

Proposition 2.1. For every cubical d-polytope {d > 3), the following are equivalent: 

• All dual manifolds of P are orientable. 

• The 2-skeleton of P has a consistent edge orientation. 

• The 2-skeleton of P contains no immersion of a cubical Mobius strip with parallel inner edges. 
2.5 From PL immersions to cubical PL spheres 

The emphasis in this paper is on cubical convex d-polytopes. In the more general setting of cu- 
bical PL (d — l)-spheres, one has more flexible tools available. In this setting, Babson & Chan [3] 
proved that "all PL codimension 1 normal crossing immersions appear." The following sketch 
is meant to explain the Babson-Chan theorem geometrically (it is presented in a combinatorial 
framework and terminology in [3] ) , and to briefly indicate which parts of their construction are 
available in the polytope world. 
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Construction 1: Babson-Chan [3] 



Input: A normal crossing immersion j : A4'^ ^ — > ^ of a triangulated PL manifold 
q£ dimension d — 2 into a PL simplicial (d — l)-spliere. 

Output: A cubical PL (d — l)-spliere with a dual manifold immersion PL-equivalent to j. 

(1) Perform a barycentric subdivision on and S'^~^. 




Figure 6: Step 1. Performing a barycentric subdivision. (We illustrate the impact 
of the construction on 2-ball, which might be part of the boundary of a 2-sphere. The 
immersion which is shown in bold has a single double-intersection point.) 

(Here each i-simplex is replaced by {i + 1)! new i-simplices, which is an even num- 
ber for i > 0. This step is done only to ensure parity conditions on the /-vector, 
especially that the number of facets of the final cubical sphere is congruent to the 
Euler characteristic of Barycentric subdivisions are easily performed in the 

polytopal category as well, see Ewald & Shephard [11].) 

(2) Perform a "cubical barycentric subdivision" on A4'^~^ and S'^~^. 




Figure 7: Step 2. Performing a cubical barycentric subdivision. 

(This is the standard tool for passage from a simplicial complex to a PL-homeomorphic 
cubical complex; here every i-simplex is subdivided into i + 1 different i-cubes. Such 
cubations can be performed in the polytopal category according to Shephard [29]: 
If the starting triangulation of S'^~^ was polytopal, the resulting cubation will be 
polytopal as well.) 

(3) "Thicken" the cubical {d — l)-sphere along the immersed (d — 2)-manifold, to obtain 
the cubical (d - l)-sphere BC{S'^-^,j{M'^-^)). 

(In this step, every (d — 1 — i)-cube in the i-fold multiple point locus results in a new 
{d— l)-cube. The original immersed manifold, in its cubified subdivided version, now 
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appears as a dual manifold in the newly resulting (d — l)-cubes. This last step is the 
one that seems hard to perform for polytopes in any non-trivial instance.) 




Figure 8: The outcome of the Babson-Chan construction: A cubical sphere with a 
dual manifold immersion that is PL-equivalent to the input immersion j. 



3 Lifting polytopal subdivisions 
3.1 Regular balls 

In the following, the primary object we deal with is a regular ball: a regular polytopal subdivi- 
sion B of a convex polytope P = \B\. 

Definition 3.1 (regular subdivision, lifting function). A polytopal subdivision B is regular (also 
known as coherent or projective) if it admits a lifting function, that is, a concave function 
/ : |P| — > M whose domains of linearity are the facets of the subdivision. (A function g : D ^ M 
is concave if for all x,y e D and < A < 1 we have g{Xx + (1 — X)y) > Xg{x) + (1 — X)g{y).) 

In this definition, subdivisions of the boundary are allowed, that is, we do not necessarily require 
that the faces of P = \B\ are themselves faces in B. 

In the sequel we focus on regular cubical balls. Only in some cases we consider regular non- 
cubical balls. 

Example. If (P, F) is an almost cubical polytope, then the Schlegel diagram based on F, which 
we denote by Schlegel(P, F), is a regular cubical ball (without subdivision of the boundary). 

Lemma 3.2. If B is a regular cubical d-ball, then there is a regular cubical ball B' without 
subdivision of the boundary, combinatorially isomorphic to B. 

Proof. Using a positive lifting function / : |S| — s- M, the d-ball B may be lifted to B in M'^"'"-'^, by 
mapping each X G \B\ to {x,f{x)) G 1^^^+^ 

Viewed from p := Xe^+i for sufficiently large A, this lifted ball will appear to be strictly convex, 
that is, its boundary is a convex polytope (rather that a boundary subdivision of a convex 
polytope). Thus one may look at the polytopal complex that consists of the cones spanned 
by faces of B with apex p. This polytopal complex is regular, since it appears convex when 
viewed from p, which yields a lifting function for the restriction of B to the hyperplane given by 
Xd+i = 0, which may be taken to be B' . □ 
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p 




B 



Figure 9: Illustration of the 'convexification' of a regular ball (Lemma 3.2). 

3.2 Lifted balls 

When constructing cubical complexes we often deal with regular cubical balls which are equipped 
with a lifting function. A lifted d-hall is a pair {B, h) consisting of a regular d-ball B and a lifting 
function h of B. The lifted boundary of a lifted ball {B, h) is the pair {dB, h\QQ). 

If {B^h) is a lifted d-ball in W^' then lift(S, /i) denotes the copy of B in W^'^^ with vertices 
{v, h{v)) € V ^ vert(S). (In the sequel we sometimes do not distinguish between these two 

interpretations of a lifted ball.) We rely on Figure 10 for the illustration of this correspondence. 




Figure 10: A lifted cubical ball {B,h) and its lifted copy lift(B, ft). The figure on the right shows the convex 
hull Q = conv(lift(B, h)). 

Notation. We identify M*^ with x {0} C M'^^^, and decompose a point x G M'^+^ as a; = 
(7r(a;), 7(2;)), where 'y{x) is the last coordinate of x and tt : M.'^'^^ M"' is the projection that 
eliminates the last coordinate. 

Often a lifted ball {B, ip) is constructed as follows: Let P be a d-polytope (in R'^) and Q C W^'^^ 
a (d+ l)-polytope such that t^{Q) = P. Then the complex B' given as the set of upper faces of Q 
determines a lifted polytopal subdivision {B, tp) of P (where B := '7t{B') and ^p is determined the 
vertex heights 7(1'), v G vert(;B')). Hence lift(;B,V') equals B'. Compare again Figure 10. 

The lifted boundary subdivision of a d-polytope P is a pair (5''"^,^') consisting of a polytopal 
subdivision S'^~^ of the boundary of P and a piece-wise linear function ip : \dP\ — > M such that 
for each facet F of P the restriction of "0 to P is a lifting function of the induced subdivision 
5^^^^ n F of F. 
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3.3 The patching lemma 

Often regular cubical balls are constructed from other regular balls. The following "patching 
lemma", which appears frequently in the construction of regular subdivisions (see [22, Cor. 1.12] 
or [8, Lemma 3.2.2]) is a basic tool for this. 

Notation. For a d-polytope P C M"^', a polytopal subdivision T oi P and a hyperplane H in W^' , 
we denote hy T r\ H the restriction of T to H, which is given by 

TnH ■.= {FnH : FeT}. 

For two d-polytopes P, Q with Q C P and a polytopal subdivision T oi P we denote hy T HQ 
the restriction oi T to Q, which is given by 

TnQ:={FnQ : F G T} . 

By fac(5) we denote the set of facets of a complex S. 

Lemma 3.3 ("Patching lemma"). Let Q be a d-polytope. Assume we are given the following 
data: 

t> A regular polytopal subdivisions of Q (the "raw subdivision"). 
t> For each facet F of S, a regular polytopal subdivision Tp of F, 

such that TpHF' = Tf' nF for all facets F, F' of S . 
\> For each facet F of S, a concave lifting function hp ofTp, 

such that hp{x) = hpi{x) for all x & F D F' , where F, F' are facets of S. 

Then this uniquely determines a regular polytopal subdivision lA = yjpFp of Q (the "fine sub- 
division"). Furthermore, for every lifting function g of S there exists a small > such that 
for all e in the range Eq > e > the function g + eh is a lifting function oflA, where h is the 
piecewise linear function h : \Q\ — > R which on each F € S is given by hp. 



Proof. Let g' be a lifting function of S. For a parameter e > we define a piece-wise linear 
function cj);; : |P| — > M that on a; € F G fac(5) takes the value 0e(a;) = g{x) + ehp{x). (It is 
well-defined since the hp coincide on the ridges otS.) The domains of linearity of (j)^ are given by 
the facets of the "fine" subdivision lA. If e tends to zero then tends to the concave function g. 
This implies that there exists a small eq > Q such that (j)e is concave and thus a lifting function 
of U, for eo > e > 0. □ 



3.4 Products and prisms 

Lemma 3.4 ("Product lemma"). Let {Bi, hi) be a lifted cubical di-ball in M'^'i and {B2, h2) be a 
lifted cubical d2-ball in 

Then the product Bi x B2 of Bi and B2 is a regular cubical {di + d2)-ball in M'^i+'^a. 

Proof. Each cell of Bi x is a product of two cubes. Hence Bi x is a cubical complex. A lifting 
function h oi Bi x B2 is given by the sum of hi and h2, that is, by h{{x,y)) := hi{x) + h2{y), 
for X G \Bi\,y € \B2\. □ 

As a consequence, the prism prism(C) over a cubical d-complex C yields a cubical (d + 1)- 
dimensional complex. Furthermore, the prism over a regular cubical ball B yields a regular 
cubical {d + l)-ball. 
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3.5 Piles of cubes 



For integers £i, . . . , £rf > 1, the pile of cubes Pd{h, • • • , ^d) is the cubical d-bah formed by ah unit 
cubes with integer vertices in the d-polytope P := x . . . x [0, ^rf], that is, the cubical d-ball 

formed by the set of all d-cubes 

C{ki, . . . ,kd) := [ki,ki + 1] X ■ ■ ■ X [kd, kd + 1] 

for integers < ki < £i together with their faces [33, Sect. 5.1]. 

The pile of cubes Pdi^i, ■ ■ ■ Ad) is a product of 1-dimensional subdivisions, which are regular. 
Hence the product lemma implies that Pd{(-i, ■ ■ ■ i^d) is a regular cubical subdivision of the 
d-polytope P. 

3.6 Connector polytope 

The following construction yields a "connector" polytope that may be used to attach cubical 
4-polytopes resp. regular cubical 4-balls without the requirement that the attaching facets are 
projectively equivalent. 

Lemma 3.5. For any combinatorial 3-cube F there is a combinatorial A-cube C that has both (a 
projective copy of) F and a regular 3-cube F' as (adjacent) facets. 

Proof. After a suitable projective transformation we may assume that F C has a unit 
square Q as a face. Now the prism F x I over F has F and Q x I as adjacent facets, where the 
latter is a unit cube. □ 

4 Basic construction techniques 
4.1 Lifted prisms 

While there appears to be no simple construction that would produce a cubical (d + l)-polytope 
from a given cubical d-polytope, we do have a simple prism construction that produces regular 
cubical {d + l)-balls from regular cubical d-balls. 

Construction 2: Lifted prism 



Input: A lifted cubical d-ball {B, h). 

Output: A lifted cubical {d + l)-ball LiftedPrism(;B, h) 

which is combinatorially isomorphic to the prism over B. 

We may assume that the convex lifting function h defined on P := \B\ is strictly positive. 
Then the hfted facets of LiFTEDPRiSM(;i3, h) may be taken to be the sets 

F := {{x,t,h{x)):xeF,-h{x)<t<+h{x)}, Fefac{B). 



If B does not subdivide the boundary of P, then LiftedPrism(;B, h) does not subdivide the 
boundary of | LiftedPrism(;B, h)\. In this case P := | LiftedPrism(S, h)\ is a cubical {d+ 1)- 
polytope whose boundary complex is combinatorially isomorphic to the boundary of the prism 
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over B. The /-vector of P is then given by 




LiftedPrism(;B, h) 



Figure 11: The lifted prism of a lifted cubical d-ball {B, h), displayed for d = 2. The result is a (regular) cubical 
(d + l)-ball that is combinatorially isomorphic to the prism over B. 

Proposition 4.1 (Dual manifolds). Up to PL-homeomorphism, the cubical ball LiftedPrism(B, h) 
has the following dual manifolds: 

• J\f X I for each dual manifold M of B, 

• one {d — 1)- sphere combinatorially isomorphic to dB. □ 
4.2 Lifted prisms over two balls 

Another modification of this construction is to take two different hfted cubical balls {Bi , hi ) and 
{B2,h2) with the same lifted boundary complex (that is, dBi = dB2 with hi{x) = h2{x) for all 
X G dBi = dB2) as input. In this case the outcome is a cubical {d+ l)-polytope which may not 
even have a cubification. 

Construction 3: Lifted prism over two balls 



Input: Two lifted cubical d-balls and (S2,/i2) 

with the same lifted boundary. 

Output: A cubical {d + l)-polytope LiftedPrism((Si, /ii), /12)) 
with lifted copies of Bi and Bi in its boundary. 

If both balls do not subdivide their boundaries, we set B'^ := B^ and /i'^ := h^ for k G {1, 2}. 
Otherwise we apply the construction of the proof of Lemma 3.2 simultaneously to both lifted 
cubical balls and {B2,h2) to obtain two lifted cubical d-balls {B'i,h'i) and (-i32,/i2) 

with the same support Q = = |,B2| which do not subdivide the boundary of Q. 

We can assume that h[, h'2 are strictly positive. Then Q := LiftedPrism((;Bi, /ii), {B2, /i2)) 
is defined as the convex hull of the points in 

{{x,+h'i{x)):x(;i\B[\} U {{x,-h'2{x)):x(^\B'2\}, 
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Since B'^ and 82 both do not subdivide their boundaries, each of their proper faces yields a 
face of Q. Furthermore, Q is a cubical {d + l)-polytope whose /-vector is given by 



= f/o(Si)+/o(^2) forA: = 0, 

I fkiBi) + fk{B2) + fk-i{dBi) ioiO <k<d. 



See Figure 12. 




LiftedPrism((^i, hi), {B2, 

Figure 12: The lifted prism over two lifted cubical d-balls {Bi, hi) and {B2, /12), displayed for d — 2. The outcome 
is a cubical (d + l)-polytope. 

4.3 Schlegel caps 

The following is a projective variant of the prism construction, applied to a d-polytope P. 
Construction 4: Schlegel cap 

Input: An almost cubical d-polytope {P,Fq) 

Output: A regular cubical d-bah SchlegelCap(P, Fq), with P C | SchlegelCap(P, _Fo)| 
which is combinatorially isomorphic to the prism over Schlegel(P, i^). 

The construction of the Schlegel cap depends on two further pieces of input data, namely 
on a point xq G beyond Fq (and beneath all other facets of P; cf. [15, Sect. 5.2]) and 
on a hyperplane H that separates xq from P. In terms of projective transformations it is 
obtained as follows: 

(1) Apply a projective transformation that moves Xq to infinity while fixing H point- 
wise. This transformation moves the Schlegel complex C{dP)\{FQ} to a new cubical 
complex £. 

(2) Reflect the image £ of the Schlegel complex in H, and call its reflected copy £' . 

(3) Build the polytope bounded by £ and £' . 

(4) Reverse the projective transformation of (1). 



12 



p 

H 



H 



Figure 13: Construction steps of the Schlegel cap over an almost cubical polytope. 

An alternative description, avoiding projective transformations, is as follows: 

(1) For each point x in the Schlegel complex C{dP)\{FQ\ let x be the intersection point 
of H and the segment [xq^ x\, and let a;' be the point on the segment [xq, x] such that 

] form a harmonic quadruple (cross ratio —1). 
That is, if 17 is a direction vector such that x = Xq + tv for some t > 1 denotes the 
difference x — xq, while x = xq + v lies on H, then x' = xq+ 2^^v. 

(2) For each face G of the Schlegel complex, G' := {x' : x G G} is the "projectively 
reflected" copy of G on the other side of H. 

(3) The Schlegel cap SchlegelCap(P, Fq) is the regular polytopal bah with faces G, G' 
and conv(G U G') for faces G in the Schlegel complex. 




Figure 14: Constructing the Schlegel cap via cross ratios. 



5 A small cubical 4-polytope with a dual Klein bottle 

In this section we present the first instance of a cubical 4-polytope with a non-orientable dual 
manifold. By Proposition 2.1 this instance is not edge-orientable. Hence, its existence also 
confirms the conjecture of Hetyei [17, Conj. 2, p. 325]. Apparently this is the first example of a 
cubical polytope with a non-orientable dual manifold. 

Theorem 5.1. There is a cubical 4-polytope P72 with f -vector 

fiPn) = (72, 196, 186, 62), 
one of whose dual manifolds is an immersed Klein bottle of f -vector (80, 160, 80). 



13 



Step 1. We start with a cubical octahedron Og, the smaUest cubical 3-polytope that is not a 
cube, with /-vector 

/(Og) = (10,16,8). 







'/ / 








/ 1 


1 





Figure 15: The cubical octahedron Og positioned in with a regular square base facet Q and acute dihedral 
angles at this square base. A part of one dual manifold is highlighted. 

We may assume that Og is already positioned in with a regular square base facet Q and acute 
dihedral angles at this square base; compare the figure below. The /-vector of any Schlegel 
diagram of Og is 

/(SCHLEGEL(Og,Q)) = (10,16,7). 

Let Og be a congruent copy of Og, obtained by reflection of Og in its square base followed by a 
90° rotation around the axis orthogonal to the base; compare the figure below. This results in 
a regular 3-ball with cubical 2-skeleton. Its /-vector is 

/(S2) = (16, 28, 15, 2). 

The special feature of this complex is that it contains a cubical Mobius strip with parallel inner 
edges of length 9 in its 2-skeleton, as is illustrated in the figure. 




Figure 16: The outcome of step 1 of the construction: The 2-cubical convex 3-ball Bi which contains a Mobius 
strip with parallel inner edges in the 2-skeleton. 
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Step 2. Now we perform a Schlegel cap construction on Og, based on the (unique) facet F 
of Os that is not contained in the Mobius strip mentioned above, and that is not adjacent to 
the square glueing facet Q. This Schlegel cap has the /-vector 

fiSr) = (20, 42, 30, 7), 

while its boundary has the /-vector 

fidSr) = (20, 36, 18). 

Step 3. The same Schlegel cap operation may be performed on the second copy Og. Joining 
the two copies of the Schlegel cap results in a regular cubical 3-ball B14 with /-vector 

fiBu) = (36, 80, 59, 14) 

whose boundary has the /-vector 

fidBu) = (36, 68, 34). 

The ball 814^ again contains the cubical Mobius strip with parallel inner edges of length 9 as an 
embedded subcomplex in its 2-skeleton. Compare Figure 17. 




Figure 17: The outcome of step 2 of the construction: The cubical convex 3-ball B14 which contains a Mobius 
strip with parallel inner edges in the 2-skeleton. 

Step 4. Now we build the prism over this regular cubical ball, resulting in a regular cubical 
4-ball B whose /-vector is 

f{B) = (72, 196, 198, 87, 14) 

and whose support is a cubical 4-polytope P72 := l-BI with two copies of the cubical Mobius strip 
in its 2-skeleton. Its /-vector is 

/(P72) = (72, 196, 186, 62). 

A further (computer-supported) analysis of the dual manifolds shows that there are six dual 
manifolds in total: one Klein bottle of /-vector (80, 160, 80), and five 2-spheres (four with /- 
vector (20, 36, 18), one with /-vector (36, 68, 34)). All the spheres are embedded, while the 
Klein bottle is immersed with five double-intersection curves (embedded 1-spheres) , but with no 
triple points. □ 
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6 Constructing cubifications 



A lot of construction techniques for cubifications (see Section 2.2) are available in the CW 
category. In particular, every cubical CW {d — l)-sphere S'^~^ with an even number of facets 
admits a CW cubification, that is, a cubical CW ci-ball with boundary S'^~^, according to 
Thurston [31], Mitchell [23], and Eppstein [10]. 



6.1 The Hexhoop template 



Yamakawa &; Shimada [32] have introduced an interesting polytopal construction in dimension 
3 called the Hexhoop template; see Figure 18. 




Figure 18: The Hexhoop template of Yamakawa & Shimada [32]. 



Their construction takes as input a 3-polytope P that is affinely isomorphic to a regular 3-cube, 
a hyperplane H and a cubical subdivision S of the boundary complex of P such that S is 
symmetric with respect to H and H intersects no facet of S in its relative interior. For such a 
cubical PL 2-sphere S the Hexhoop template produces a cubification. A 2-dimensional version 
is shown in Figure 19. 




Figure 19: A two-dimensional version of the Hexhoop template. 



6.2 The generalized regular Hexhoop — overview 

In the following we present a generalized regular Hexhoop construction. It is a generalization of 
the Hexhoop template in several directions: Our approach admits arbitrary geometries, works 
in any dimension, and yields regular cubifications with "prescribed heights on the boundary" 
(with a symmetry requirement and with the requirement that the intersection of the symmetry 
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hyperplane and the boundary subdivision is a subcomplex of the boundary subdivision). 
Figure 20 displays a 2-dimensional cubification (of a boundary subdivision 5 of a 2-polytope 
such that S is symmetric with respect to a hyperplane H) obtained by our construction. 




Figure 20: A cubification of a boundary subdivision of a pentagon, produced by our generalized regular Hexhoop 
construction. 

Not only do we get a cubification, but we may also derive a symmetric lifting function for 
the cubification that may be quite arbitrarily prescribed on the boundary. The input of our 
construction is a lifted cubical boundary subdivision {S'^~^,Tp) of a (i-polytope P, such that 
both P and are symmetric with respect to a hyperplane H. 




H 



Figure 21: An input for tlie generalized regular Hexhoop construction. 

Our approach goes roughly as follows. 

(1) We first produce a (d + l)-polytope T that is a symmetric tent (defined in Section 6.3) 
over the given lifted boundary subdivision (5, ^/i) of the input d-polytope P. Such a tent 
is the convex hull of all 'lifted vertices' {v,ip{v)) € R'^"'"^, v G vert(5), and of two apex 
points p^,pp>; compare Figure 22. 



Pl PR 




H 



Figure 22: A symmetric tent over tfie lifted boundary subdivision {S,'ip) of the input d-polytope P. 
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(2) Truncate T by a hyperplane H' parallel to aff(P) = M*^ c M'' x {0} that separates the 
lifted points from the apex points, and remove the upper part. 

(3) Add the polytope R := cone(p^, Q) fl cone(p^, Q) n H'j^^ where H'j^ is the halfspace with 
respect to H that contains and p^. Compare Figure 23. 

(4) Project the upper boundary complex of the resulting polytope to W^. 

The figures in this section illustrate the generalized regular Hexhoop construction for the 2- 
dimensional input polytope of Figure 21; the generalized Hexhoop construction for d = 2 yields 
2-dimensional complexes in M'^. The extension to higher dimensions is immediate, and the case 
d = 3 is crucial for us (see Section 7). It is, however, also harder to visualize: A 3-dimensional 
generalized regular Hexhoop cubification is shown in Figure 29. 




Figure 23: Sketch of the generalized regular Hexhoop construction. 

6.3 Symmetric tent over a lifted boundary subdivision 

Let P be a d-polytope that is symmetric with respect to a hyperplane H in W^. Choose a 
positive halfspace Hj^ with respect to H. Let (5, "0) be a lifted boundary subdivision of P such 
that S r\ H \s a, subcomplex of S. Define H := H -\- Me^+i, which is a symmetry hyperplane 
for P C E'^+^ The positive halfspace oi H d wiU be denoted by H+. 

The symmetric tent over (5, ip) is the lifted polytopal subdivision (T, (j)) of P given by the upper 
faces of the polytope 

T := conv(PU{Pi,p^}) 

ifpi,Pf( G M.'^'^^ are two apex points in R'^"'"^ that are symmetric with respect to the hyperplane 
H, and the upper facets of T are 



18 



• pyramids with apex point over facets F of lift(5, V') such that vr(F) C 

• pyramids with apex point over facets F of hft(5, ij:) such that ■k{F) C and 

• 2-fold pyramids with apex points p^,pj^ over ridges R of hft(5, "0) with it{R) C 

(This requires that aff (P) and vr(p^) G rehnt(P n H^).) 
Lemma 6.1. Assume we are given the following input. 

P a convex d-polytope in M'^, 

{S,ip) a lifted boundary subdivision of P, 

H a hyperplane in such that 

• P and {S, ip) are both symmetric with respect to H , and 

• S r\H is a subcomplex of S, and 
9L' ^"^0 points in P cW^ such that 

• £ reUnt(P n and 

• QliQr ^''^^ symmetric with respect to H. 

Then for every sufficiently large height h > the {d+ l)-polytope T := conv{hft(iS, 'i/'),PL,P_R} 
with Pi := (qi, h) G Hj^ and p^ := (q^, h) ^ Hj^ is a symmetric tent over (5, ■0)- 

This can be shown for instance by using the Patching Lemma (Lemma 3.3). 



6.4 The generalized regular Hexhoop in detail 

In this section we specify our generahzation of the Hexhoop template and prove the following 
existence statement for cubifications. 

Theorem 6.2. Assume we are given the following input. 
P a convex d-polytope in , 

{S'^~^,'il)) a lifted cubical boundary subdivision of P, and 
H a hyperplane in M'^ such that 

• P and {S'^~^,'ip) are symmetric with respect to H, and 

• S'^~^ n H is a subcomplex of S'^~^ . 

Then there is a lifted cubification {B'^,(t>) of {S'^~^,'ip). 
The proof relies on the following construction. 

Construction 5: Generalized regular Hexhoop 



Input: 

P a convex d-polytope P in W^. 

{S'^"^, if:) a lifted cubical boundary subdivision of P. 

H a hyperplane in such that 

• P and {S'^~^,il)) are symmetric with respect to H, and 

• S'^~^ n i/ is a subcomplex of S'^~^ . 

Output: 

{B'^,(j)) a symmetric lifted cubification of (5"^""^,^) given by a cubical d-ball C 
in M'^+i. 



(1) Choose a positive halfspace i7+ with respect to H, and a point € relint(P n 

Define q^j := pf^^ , where the upper index denotes the mirrored copy with respect 
to H = H + Red+i- 
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By Lemma 6.1 there is a height h > such that 

T:=conv{m{S''~\ij),PL,PR} 

with := (qi, h) and p^ := {q^, h) forms a symmetric tent over {S'^~^ , if)) . 

(2) Choose a hyperplane H' parahel to aff (P) C that separates {Pl^Pr} ^iid hft(5'^~^, tp). 
Let H'j^ be the halfspace with respect to H' that contains p^ and p^. 




H 



Figure 24: Step 2. The hyperplane H' separates {p^jP^} from lift(5'' 

(3) Define the "lower half" of the tent T as 

r_ := TnHL, 
whose "top facet" is the convex d-polytope Q := T (1 H'. 




Figure 25: Step 3. The "lower half" T- of T. 

(4) Define the two d-polytopes 

Ql ■= conv{?; G vert((5) : v G -f^+}, 
Qr := conv {i; G vert{Q) : v G 

Let Fl := H' n conv(p^, P n H), the unique facet of that is not a facet of Q. 




Figure 26: Step 4. Define and Qa. 
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(5) Construct the polytope 

R := cone{p^,Q) fl cone(p^, Q) fl H'^ 



H = H + Mcd+i 




Figure 27: Step 5. The polytope 7? := cone(pj^, Q) D cone(p^, Q) n H'^. 

The complex C in question is given by the upper facets of the {d + l)-polytope 

U := T^UR. 



Lemma 6.3 (Combinatorial structure of Q). 

The vertex set of Q consists of 

• the points conv(p£^, v) R H' for vertices v € vert(hft(5, such that ■7r{v) C H^, and 

• the points conv(p^,t;) n H' for vertices v G vert(hft(5, ■0)) such that tt{v) C H^. 

The facets of Q are 

(a) the combinatorial cubes conv{p]^, F) n H' for facets F ofli^t{S,^p) such that F C -ff+, 

(b) the combinatorial cubes conv(p^,F) fl H' for facets F of lih{S , ip) such that F C 

(c) the combinatorial cubes conv(p^,pj:j, F) n H' for (d — 2)-faces F of hft(5, ijj) with F C H. 

Proof. By the definition of a symmetric tent, upper facets of the symmetric tent T are 

• the pyramids with apex point Pi over facets F of hft(iS, ^p) such that F C 

• the pyramids with apex point pj^ over facets F of hft(5, ■0) such that F C and 

• the 2-fold pyramids with apex points Pi^,pji over ridges R of lift(5,'0) with R C H. 

Since Q is the intersection of T with H, the polytope Q has the vertices and facets listed above. 
It remains to show that the facets of type (c) are combinatorial cubes. Let F he a {d — 2)-face 
of lift(5, ip) such that F C H. Every point on the facet lies in the convex hull of F with a 
unique point on the segment [p^ , p^] . Thus the facet is combinatorially isomorphic to a prism 
over F. □ 
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Let a d- dimensional half-cube be the product of a combinatorial {d — 2)-cube and a triangle. A 
combinatorial half-cube is a polytope combinatorially isomorphic to a half-cube. 



Lemma 6.4 (Combinatorial structure of T_). 

The vertices of T- are the vertices o/lift(5, V') and the vertices of Q. Furthermore, the upper 
facets ofT^ are 

(a) the combinatorial cubes cone{pi, F) riH'_ n (M'^ x R+) for facets F of Q such that F C H^, 

(b) the combinatorial cubes cone(p^, F) r\H'_ H {W^ x for facets F of Q such that F C H_, 

(c) the combinatorial half-cubes cone{pi, F) Ci cone{pji, F) H H'_ for facets R of Q that inter- 
sect H , and 

(d) Q. 

The facet defining hyperplanes of the upper facets of T_ are 

(a) aff(p^,F) for facets F of Q such that F C H^, 

(b) a.S{pji,F) for facets F of Q such that F C 

(c) aff(p£,p^,F) for facets F of Q that intersect H, and 

(d) afr(Q). 

Proof. Since T_ is the intersection of T with H'_, the upper facets of T_ are given by Q plus 
the intersections of the upper facets of T with H'_, and the vertices of r_ are the vertices of T 
and the vertices of Q. □ 

Lemma 6.5 (Combinatorial structure of R). 

The set of vertices of R consists of the vertices of Q and all points in V" := vert(i?) \ vert(Q). 
Furthermore, the set of (all) facets of R consists of 

(a) the combinatorial cubes conv(p^,F) n for facets F of Q such that F C H^, 

(b) the combinatorial cubes conv(p£,F) CiH^ for facets F of Q such that F C 

(c) the combinatorial half-cubes conv(p^,F) n conY{p]^, F) for facets F of Q that intersect H, 
and 

(d) Q. 

The set of facet defining hyperplanes of the facets of R consists of 

(a) aff(p^,F) for facets F of Q such that F C 

(b) aS.{pi,F) for facets F of Q such that F C ^ 

(c) a&{pi^,pji, F) for facets F of Q such that F intersects H, and 

(d) aff(Q). 

Proof of Theorem 6. 2. We show that the complex C' given by the upper facets of the polytope U 
of Construction 5 determines a lifted cubification {B'^, cj)) of {S'^~^,il)). 

First observe that no vertex of T_ is beyond a facet of R, and no vertex of R is beyond a facet 
of T_. Hence the boundary oiU = conv(T_ U R) is the union of the two boundaries of the two 
polytopes, excluding the relative interior of Q. 

Define the vertex sets V := vert(lift(5t, ^j)), V := vevt{Q) and V" := vert(fi) \ V . Then 

• each vertex of V is beneath each facet of R that is of type (a) or (b), and 

• each vertex of V" is beneath each facet of r_ that is of type (a) or (b). 

Hence these four types of facets are facets of U that are combinatorial cubes, and the set of 
vertices of U is given by the union of V,V' and V" . It remains to show that each hyperplane 
aff(p^,p^, F), where F is a facet of Q that intersects H, is the affine hull of a cubical facet 
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of U . To see this, observe that there are two facets of i?, T_ respectively, that are both 

contained in the affine huh of F. These two facets F- are both half-cubes that intersect 
in a common [d — l)-cube, namely F. Furthermore, all vertices of i<+ and of F_ that are not 
contained in aff (F) are contained in H. Hence the union of and F- is a combinatorial cube. 

Thus every upper facet of [/ is a combinatorial cube. Furthermore, ■k{R) = vr((5) and vr(T_) = 
|P|, so the upper facets of U determine a lifted cubical subdivision of {S'^~^,ip). □ 

Proposition 6.6 (Dual manifolds). Up to PL-homeomorphism, the generalized regular Hexhoop 
cubification B'^ of S'^~^ has the following dual manifolds: 

• M X I for each dual manifold J\f (with or without boundary) of Sl = S'^~^ n H^, 

• two [d — l)-spheres "around" Q, Q^'^ , respectively, where the upper index denotes the 
mirrored copy. 

Proof. The "main part" of the complex B'^ may be viewed as a prism of height 4, whose dual 
manifolds are of the form A/" x /, as well as four (d — l)-balls. This prism is then modified by 
glueing a full torus (product of the {d — 2)-sphere 5*^"^ fl TC with a square into its "waist." 
This extends the dual manifolds Af x I without changing the PL-homeomorphism type, while 
closing the four (d — l)-balls into two intersecting, embedded spheres. 

We refer to Figure 28 (case d = 2) and Figure 29 {d = 3) for geometric intuition. 




(a) Dual product manifolds N x I (b) Dual spheres 

Figure 28: The dual manifolds of a 2-dimensional generalized regular Hexhoop. 




Figure 29: A 3-dimensional cubification produced by the generalized regular Hexhoop construction. For every 
embedded dual circle N which intersects -ff+ \-ff and H- \H , there is an embedded dual 2-ball M' with boundary N 
in the the cubification. (This is a cubification for the case "singleS" introduced in Section 7.) 

□ 
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7 Cubical 4-polytopes with prescribed dual manifold immersions 



Now we use our arsenal of cubical construction techniques for the construction of cubifications 
with prescribed dual 2-manifold immersions, and thus approach our main theorem. 

For this we ask for our input to be given by normal crossing PL-immersions whose local geometric 
structure is rather special: We assume that A4'^~^ is a (d— l)-dimensional cubical PL-manifold, 
and J : TW^""*^ S-^ M'^ is a grid immersion, a cubical normal crossing codimension one immersion 
into M'^ equipped with the standard unit cube structure. 

7.1 From PL immersions to grid immersions 

In view of triangulation and approximation methods available in PL and differential topology, 
the above assumptions are not so restrictive. (See, however, Dolbilin et al. [9] for extra problems 
and obstructions that may arise without the PL assumption, and if we do not admit subdivisions, 
even for the high codimension embeddings/immersions.) 

Proposition 7.1. Every locally flat normal crossing immersion of a compact (d — l)-manifold 
into is PL-equivalent to a grid immersion of a cubification of the manifold into the standard 
cube subdivision ofW^. 

Proof, may replace any PL-immersion of A^'^"^ by a simplicial immersion into a suitable trian- 
gulation of M'^. The vertices of may be perturbed into general position. 

Now we overlay the polyhedron j{A4'^~^) with a cube structure of M*^ of edge length e for suitably 
small e > 0, such that the vertices of j{A4'^~^) are contained in the interiors of distinct d-cubes. 

Then working by induction on the skeleton, within each face of the cube structure, the restriction 
of to a k-face — which by local flatness consists of one or several {k — l)-cells that 

intersect transversally — is replaced by a standard cubical lattice version that is supposed to 
run through the interior of the respective cell, staying away distance e' from the boundary of the 
cell; here we take different values for e' in the situation where the immersion is not embedded 
at the vertex in question, that is, comes from several disjoint neighborhoods in Jv['^~^ . 




Figure 30: Illustration of the proof of Proposition 7.1. 



The resulting modified immersion into will be cellular with respect to a standard cube 
subdivision of edge length jtE for a suitable large A^. Figure 30 illustrates this for d = 2. □ 
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7.2 Vertex stars of grid immersions of surfaces 

From now on, we restrict our attention to the case of d = 3, that is, 2-manifolds and 4-polytopes. 

There are nine types of vertex stars of grid immersions of surfaces, namely the following five 
vertex stars of a regular vertex, 



3 a a --m- 




singles single4a single4b singleS single6a single6b 

plus two vertex stars with double intersection and the vertex star of a triple intersection point: 

-e- -s- ifr 

doubleSa doubleSb triplel2 

For the constructions below we will require that the grid immersion j : Ai"^ S-^ that we start 
out with is locally symmetric, that is, that at each vertex w of j(A^^) there is plane H through 
w such that for each vertex v with j{v) = w the image of the vertex star of v is symmetric 
with respect to H. Thus we require that H is a symmetry plane separately for each of the (up 
to three) local sheets that intersect at w. Such a plane H is necessarily of the form Xi = k, 
or Xi — Xj = k. In the first case we say H is a. coordinate hyperplane, and in other 
cases it is diagonal. 

Proposition 7.2. Any grid immersion of a compact cubical 2-manifold into M'^ is equivalent to 
a locally symmetric immersion of the same type. 

Proof. All the vertex stars displayed above satisfy the local symmetry condition, with a single 
exception, namely the star "single6b" of a regular vertex with six adjacent quadrangles. As 
indicated in Figure 31, a local modification of the surface solves the problem (with a suitable 
refinement of the standard cube subdivision) . □ 




single6b 

Figure 31: Local modification used to "repair" the case "singleGa." 
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7.3 Main theorem (2-manifolds into cubical 4-polytopes) 

Theorem 7.3. Let j: S-> be a locally flat normal crossing immersion of a compact 2- 
manifold (without boundary) A4 into M"^. 

Then there is a cubical A-polytope P with a dual manifold M.' and associated immersion y : 
M' S-^ \dP\ such that the following conditions are satisfied: 

(i) M' is a cubical subdivision of M, and the immersions j (interpreted as a map to M.^ U 
{oo} = S^) and y are PL-equivalent. 

(ii) The number of facets of P is congruent modulo 2 to the number t{j) of triple points of the 
immersion j. 

(iii) // the given surface Ai is non-orientable and of odd genus, then the cubical 4-polytope P 
has an odd number of facets. 

The core of our proof is the fohowing construction of cubical 3-balls with a prescribed dual 
manifold immersions. 



Construction 6: Regular cubical 3-ball with a prescribed dual manifold 

Input: A 2-dimensional closed (that is, compact and without boundary) cubical PL- 
surface A4 , and a locally symmetric codimension one grid immersion j : A4 

iPsih, i2,i3)\ ewe. 

Output: A regular convex 3-ball B with a dual manifold Ai' and associated immersion 
y : A4' S-^ \B\ such that the following conditions are satisfied: 

(i) Ai' is a cubical subdivision of Ai, and the immersions j and y are PL- 
equivalent. 

(ii) The number of facets of B is congruent modulo two to the number t{j) of 
triple points of the immersion j. 

(1) Raw complex. Let ^ be a copy of the pile of cubes ^3(^1 + l,i2 + I, is + 1) with all 
vertex coordinates shifted by —^1. (Hence Xj G {— ^, |, |, . . . , £j -|- ^} for each vertex 
X € vert (A).) 

Due to the local symmetry of the immersion, and the choice of the vertex coordinates 
of A, the following holds: 

i> Each vertex of j{Al) is the barycenter of a 3-cube C of A. 

> For each 3-cube C of ^ the restriction {C,j{Ai) fl C) is locally symmetric. 

(2) Local subdivisions. We construct the lifted cubical subdivision i3 of ^ by induction 
over the skeleton: For k = 1,2,3, B^ will be a lifted cubical subdivision of the A:- 
skeleton J^k{A), with the final result B := B^. For each A:-face G .A we take for the 
restriction B^ H F a congruent copy from a finite list of templates. 

Consider the following invariants (for k € {1, 2, 3}). 

(Ifcl) Consistency requirement. 

For every A;-face Q € Tk{A) and every facet F of Q, the induced subdivision 
B'TiF equals B^''^ n F. 
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(Ifc2) PL equivalence requirement. 

For every A;- face Q G ^jt(^) and every dual manifold J\f of Q (with bound- 
ary) the cubical subdivision B^nQ has a dual manifold that is PL-equivalent 
to j{Af) n Q. 

(Ifc3) Symmetry requirement. 

Every symmetry of {Q, j{Ai)nQ) for a A:-face Q G ^fc(^) that is a symmetry 
of each sheet of j{A4) n Q separately is a symmetry of {Q , B'' H Q) . 

(Ifc4) Subcomplex requirement. 

For every diagonal symmetry hyperplane Hq of a facet Q of ^ and every 
facet F of Q the (lifted) induced subdivision n {F Ci H) is a (lifted) 
subcomplex of S'^. 

These invariants are maintained while iteratively constructing B^ and B^. The result- 
ing lifted cubical subdivision B"^ of A will satisfy (I3I) and (I32), but not in general 
the other two conditions. 

Subdivision of edges. Let e be an edge of A. 

• If e is not intersected by the immersed manifold, then we subdivide the edge by an 
affine copy Bl of the following lifted subdivision U2 ■= {U2, h) of Pi{2): 



• • • 

1 

4 



If e is intersected by the immersed manifold, then we subdivide the edge by an 
affine copy Bl of the following lifted subdivision U3 := {U^, h) of -Pi (3): 



1 1 

6 6 



Observe that (Iil)-(Ii4) are satisfied. 



Subdivision of 2-faces. Let Q be a quadrangle of A, and w the unique vertex of j(-^) 
that is contained in Q. There are four possible types of restrictions of the grid immer- 
sion to Q: 




empty single2a single2b double 



(a) In the cases "single2a" and "double" there is a coordinate hyperplane H such that 
{Q.,j{M) n Q) is symmetric with respect to H, and a vertex v of M. such that 
j[v) = w and the image of the vertex star is contained in H. Let F be a facet 
of Q that does not intersect H. Then Bq is taken to be a copy of the product 
{B^ n F) X U3: 
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(b) If the immersion does not intersect Q, then B'q is a copy of the hfted cubical 
2-complex V which arises as the cubical barycentric subdivision of the stellar 
subdivision of [— |, ^]^: 





(c) In the case "single2b" we define 0q as an affine copy of the lifted cubical 2- 
complex V', which is given by V truncated by four additional planes: 





Observe that the conditions (I2l)-(l24) are satisfied. 
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(5) Subdivision of 3-cubes. Let Q be a facet of A and w the unique vertex of j{M) 
that is mapped to the barycenter of Q. Let S := B'^ H Q he the induced hfted cubical 
boundary subdivision of Q. 

All templates for the lifted cubification Bq of S arise either as a generalized regular 
Hexhoop, or as a product of with a lifted cubical subdivision of a facet of Q. 

(a) For the following four types of vertex stars we use a product with U-^: 

% ^ m 

single4a doubleSa doubleSb triple 12 

In all these cases there is a coordinate symmetry plane H such that H Q is 
a sheet of j{M.) H Q. Hence all facets of Q that intersect H are subdivided by 
X or W3 X U2- Let F be one of the two facets of Q that do not intersect H. 
Then the product {B'^ DF) xUs yields the lifted subdivision Bq of Q. Clearly Bq 
is consistent with (I3I) and (I32). 

(b) In the remaining five cases we take a generalized regular Hexhoop with a diagonal 
plane of symmetry of Q to produce Bq: 



□1 a 




empty singleS single4b singleS singlc6a 

In each of these cases, {Q,Q r)j{A4)) has a diagonal plane H of symmetry. This 
plane intersects the relative interior of two facets of Q. Since (l24) holds, no facet 
of iS = n Q intersects H in its relative interior. By (l23) the lifted boundary 
subdivision S is symmetric with respect to H. Hence all preconditions of the 
generalized regular Hexhoop construction are satisfied. The resulting cubification 
Bq satisfies (I3I) and (I32). 



7.4 Correctness 

Proposition 7.4. Let M be a 2- dimensional closed cubical PL-surface, and j: J^? a locally 

symmetric codimension one grid immersion. 

Then the cubical 3-ball B given by Construction 6 has the following properties: 

(i) B is regular, with a lifting function ip. 

(ii) There is a dual manifold Ai' of B and associated immersion y : A4' S-> \B\ such that M' 
is a cubical subdivision of M, and the immersions j and y are PL-equivalent. 

(iii) The number of facets of B is congruent modulo two to the number t{j) of triple points of 
the immersion j. 

(iv) There is a lifted cubification C of (55, V'las) with an even number of facets. 
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Proof, (i) Regularity. By construction the lifting functions ^l^F, F G fac(^), satisfy the consis- 
tency precondition of the Patching Lemma (Lemma 3.3). Since every pile of cube is regular the 
Patching Lemma implies that B is regular, too. 

(ii) PL- equivalence of manifolds is guaranteed by Property (I32). 

(iii) Parity of the number of facets. For each 3-cube Q of A, its cubification Bg is either a 
product Bp X U3 (where Bp is a cubification of a facet F oi Q), or the outcome of a generalized 
regular Hexhoop construction. In the latter case the the number of facets of Bq is even. In 
the first case the number of facets depends on the number of 2-faces of Bp. The number of 
quadrangles of Bp is odd only in the case "double," if n F has a double intersection point. 
Hence, fsiBq) is odd if and only if the immersion j has a triple point in Q. 

(iv) Alternative cubification. Applying Construction 6 to Psiii, 12,^3) without an immersed 
manifold yields a regular cubification C of dB with the same lifting function as B on the boundary. 
Since the immersion 9-^ M'^ has no triple points the number of facets of C is even. □ 

7.5 Proof of the main theorem 

Proof of Theorem 7. 3. Let j : A4 be a locally flat normal crossing immersion of a compact 

((i—1) -manifold M into R'^. By Proposition 7.1 and Proposition 7.2 there is a cubical subdivision 
Ai' of M with a locally symmetric, codimension one grid immersion j ': M' that is PL- 

equivalent to j. 

Construct a convex cubical 3-ball B with prescribed dual manifold immersion j' as described 
above. By Proposition 7.4(i) the ball B is regular, and by Proposition 7.4(iv) there is a cubifi- 
cation C of dB with an even number of facets and the same lifting function on the boundary. 

The lifted prism over B and C (Construction 3 of Section 4.2). This yields a cubical 4-polytope P 
with 

f3{P)=f3{B)+fs{C)+f2{dB), 

whose boundary contains B and thus has a dual manifold immersion PL-equivalent to j. 

For (ii) observe that for every cubical 3-ball the number of facets of the boundary is even. Hence 
f2{dB) is even. Since the number of facets of C is even, we obtain 

/3(P) = f3iB) = tij) mod 2. 

Now consider (iii). By a famous theorem of Banchoff [4] the number of triple points of a normal 
crossing codimension one immersion of a surface has the same parity as the Euler characteristic. 
Hence, if is a non-orientable surface of odd genus the number of triple points of j is odd, 
which implies that the cubical 4-polytope P has an odd number of facets. □ 

7.6 Symmetric templates 

The three-dimension templates constructed above, which we call the standard templates, do not 
satisfy the conditions (I33) and (I34). In particular, the symmetry requirement (I33) is violated 
by the templates corresponding to the cases "empty", "single3", and "singleGa," and it is satisfied 
by all others. For example, the standard template for "singleS" is illustrated in Figure 29; it 
satisfies (I33) since there is only one diagonal symmetry hyperplane. 

For the "empty" case an alternative template may be obtained from the cubical barycentric 
subdivision. The resulting cubification satisfies both conditions (I33) and (I34), and furthermore. 
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it has less faces — 96 facets, 149 vertices — than the standard template. 

For the case "singleS" an alternative cubification, of full symmetry, can be constructed from C" 
by truncating the lifted polytope corresponding to the lifted cubical ball C" by some additional 
hyperplanes. 

For the case "single6a" we do not know how to get a cubification of full symmetry. This is the 
main obstacle for an extension of our constructions to higher dimensions (cf. Section 11). 

8 An odd cubical 4-polytope with a dual Boy's surface 

Cubical 4-polytopes with odd numbers of facets exist by our Main Theorem 7.3. In this section 
we describe the construction of a cubical 4-polytope with an odd number of facets in more detail. 
The data for the corresponding model will be submitted to the eg-models archive. 

Theorem 8.1. There is a cubical 4-polytope P^^y with f -vector 

f = (17 718,50 784,49 599,16 533) 
that has has a Boy surface as a dual manifold. 

A grid immersion of Boy's surface. 

The construction starts with a grid immersion (cf. [25]) of Boy's surface, that is, an immersion 
of the real projective plane with exactly one triple point and three double-intersection curves in 
a pattern of three loops [7] [18] [2]. This immersion j : S-^ M'^ in shown in Figure 8. 




(a) grid immersion (b) 
Figure 32: A grid immersion the Boy's surface. Each double-intersection loop is of length four. 

The 2-manifold A4 has the /-vector f{A4) = (85, 168, 84), whereas the image of the grid 
immersion has the /-vector f{j{Ai)) = (74, 156, 84). The vertex coordinates can be chosen 
such that the image j{M.) is contained in a pile of cubes ^3(4, 4, 4). 
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A cubical 3-ball with a dual Boy's surface. 

We apply Construction 6 to the grid immersion j : S-^ to obtain a cubical 3-ball with a 
dual Boy's surface, and with an odd number of facets. 

Since the image j{A4) is contained in a pile of cubes ^3(4, 4, 4), the raw complex A given by 
Construction 6 is isomorphic to ^3(5, 5, 5). Hence we have 5"^ — 74 = 51 vertices of A that are 
not vertices of j{A4). We try to give an impression of the subdivision of the 2-skeleton of A 
in Figure 33. The /-vector of is / = (4 662, 9 876, 5 340). 

The subdivision of the boundary of A consists of 150 = 6 ■ 5 • 5 copies of the two-dimensional 
"empty pattern" template. Hence the subdivision of the boundary of A (given by n |i9^|) has 
the /-vector / = (1802, 3 600, 1800). 

The refinement B of A depends on templates that are used for the 3-cubes. We use the "sym- 
metric" templates of Section 7.6. The /-vector of B is then / = (15 915, 45 080, 43 299, 14133). 
(The "standard set" of templates yields a cubical ball with 18 281 facets.) 

Figure 34 illustrates the dual Boy's surface of the cubical 3-ball B. It has the /-vector / = 
(1998, 3 994, 1997); its multiple-intersection loops have length 16. The ball B has 612 dual 
manifolds is total (339 of them without boundary). 
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Figure 34: The dual Boy's surface of /-vector / = (1 998, 3 994, 1 997) of the cubical 3-ball B. 




A cubical 4-polytope with a dual Boy's surface 

A cubification B' of dB with an even number of facets is given by subdividing each facet of 
the raw ball A with a cubification for the empty pattern. Using the symmetric cubification for 
the empty pattern yields a regular cubical 3-ball B' with 12 000 facets. The lifted prism over B 
and B' yields a cubical 4-polytope with 27933 facets. 

However, using first stellar subdivisions, then a (simplicial) cone, and then cubical subdivisions 
on 9^3(5, 5, 5), it is possible to produce a significantly smaller alternative cubification B" of dB 
with an even number of facets. Moreover, for this one can form the cone based directly on the 
boundary complex of B and thus "save the vertical part" of the prism. The resulting cubical 
4-polytope -Psoy fo = 17718 vertices and /s = 16 533 facets. A further analysis of the dual 
manifolds of -pBoy shows that there are 613 dual manifolds in total: One dual Boy's surface of 
/-vector / = (1998, 3 994, 1997), one immersed surface of genus 20 (immersed with 104 triple 
points) with /-vector (11 470, 23016, 11 508), and 611 embedded 2-spheres with various distinct 
/-vectors. 

Verification of the instances 

All the instances of the cubical 4-polytopes described above were constructed and verified as 
electronic geometry models in the polymake system by Gawrilow & JoswiG [14], which is 
designed for the construction and analysis of convex polytopes. A number of our own tools 
for handling cubical complexes are involved as well. These cover creation, verification, and 
visualization of cubical complexes (for d E {2,3}). 

The instances are available from http://www.math.tu-berlin.de/~schwartz/c4p. 

Whereas the construction of the instances involves new tools that were writted specifically for 
this purpose, the verification procedure uses only standard polymake tools. All tools used in the 
verification procedure are parts of polymake system which have been used (and thereby verified) 
by various users over the past years (using a rich variety of classes of polytopes). 

The topology of the dual manifolds of our instances was examined using all the following tools: 

• A homology calculation code based written by Heckenbach [16]. 

• The topaz module of the polymake project, which covers the construction and analysis of 
simplicial complexes. 

• Our own tool for the calculation of the Euler characteristics. 

9 Consequences 

In this section we list a few immediate corollaries and consequences of our main theorem and of 
the constructions that lead to it. The proofs are quite immediate, so we do not give extended 
explanations here, but refer to [28] for details. 

9.1 Lattice of f-vectors of cubical 4-polytopes 

Babson &: Chan [3] have obtained a characterizazion of the Z-affine span of the /-vectors of 
cubical 3-spheres: With the existence of cubical 4-polytopes with an odd number of facets this 
extends to cubical 4-polytopes. 
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Corollary 9.1. The 'E-affine span of the f -vectors (fo, fi, f2, fa) of the cubical ^-polytopes is 

characterized by 

(i) integrality (fi € Z for all i), 

(ii) the cubical Dehn-Sommerville equations /o — /i + /2 — /s = and f2 = 8/3, and 

(iii) the extra condition fo = mod 2. 

Note that this includes modular conditions such as /2 = mod 3, which are not "modulo 2." The 
main result of Babson & Chan [3] says that for cubical d-spheres and ((i+l)-polytopes, d > 2, "all 
congruence conditions are modulo 2." However, this refers only to the modular conditions which 
are not implied by integrality and the cubical Dehn-Sommerville equations. The first example of 
such a condition is, for d = 4, the congruence (iii) due to Blind & Blind [6]. 

9.2 Cubical 4-polytopes with dual manifolds of prescribed genus 

By our main theorem, from any embedding Mg —>■ we obtain a cubical 4-polytope that has 
the orientable connected 2-manifold Mg of genus g as an embedded dual manifold. Indeed, this 
may for example be derived from a grid embedding of Mg into the pile of cubes P(l, 3, 1 + 2g). 

However, cubical 4-polytopes with an orientable dual manifold of prescribed genus can much 
more efficiently, and with more control on the topological data, be produced by means of con- 
nected sums of copies of the "neighborly cubical" 4-polytope Cf with the graph of a 5-cube 
(compare Section 2.3). 

Proposition 9.2. For each g > 0, there is a cubical 4-polytope Pg with the following properties. 

(i) The polytope Pg has exactly one embedded orientable dual 2-manifold M. of genus g with 
f -vector f {M) = {I2g + 4, 2Sg + 4, I4g + 2). 

(ii) There is a facet F of P which is not intersected by the image of the dual manifold A4, 
and which is afhnely regular, that is, there is an affine transformation between F and the 
standard cube [—1,+!]^. 

(iii) AH other dual manifolds of Pg are embedded 2-spheres. 

(iv) f{Pg) = {24g + 8, UQg + 12, 138^ + 6, 465 + 2). 

Taking now the connected sum of one example of a 4-polytope with a non-orientable dual 2- 
manifold, we obtain 4-polytopes with a non-orientable dual manifold of prescribed genus. 

Corollary 9.3. For each even g > 0, there is a cubical 4-polytope that has a cubation of the 
non-orientable connected 2-manifold Mg of genus g as a dual manifold (immersed without triple 
points and with one double-intersection curve). 

For this, one can for example construct the 4-polytope associated with the grid immersion of 
the Klein bottle of /-vector / = (52, 108, 56) as depicted in Figure 36. 

Smaller cubical 4-polytopes with non-orientable cubical 4-polytopes can be produced by means 
of connected sums of the cubical 4-polytope Pq2 of Section 5 with a dual Klein bottle, and several 
copies of the neighborly cubical 4-polytope Cf. (Some "connector cubes" of Lemma 3.5 have 
to be used.) The resulting cubical 4-polytope has rather small /-vector entries, but the set of 
multiple- intersection points consists of five double-intersection curves. 

Applying the same proof as above to the grid immersion of Boy's surface of the previous section 
yields the following result. 
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Figure 36: A grid immersion of the Klein bottle with one double-intersection curve and without triple points. 

Corollary 9.4. For each odd g > 0, there is a cubical A-polytope that has a cubation of the non- 
orientable connected 2-manifold M'g of genus g as a dual manifold (immersed with one triple 
point and three double-intersection curves of length 14). 



9.3 Higher-dimensional cubical polytopes with non-orientable dual manifolds 
Corollary 9.5. For each d > 4 there are cubical d-polytopes with non-orientable dual manifolds. 

Proof. By construction, the 4-dimensional instance Pq2 of Section 5 comes with a subdivision 
into a regular cubical 4-ball. Since one of its dual manifolds is not orientable, its 2-skeleton is 
not edge orientable, i. e. it contains a cubical Mobius strip with parallel inner edges. So if we 
now iterate the lifted prism construction of Section 4.1, then the resulting cubical d-polytopes 
{d > 4) will contain the 2-skeleton of Pq2- By Proposition 2.1 they must also have non-orientable 
dual manifolds. □ 



10 Applications to hexa meshing 

In the context of computer aided design (CAD) the surface of a workpiece (for instance a part 
of a car, ship or plane) is often modeled by a surface mesh. In order to analyze physical and 
technical properties of the mesh (and of the workpiece), finite element methods (FEM) are 
widely used. 

Such a surface mesh is either a topological mesh, that is, 2-dimensional regular CW complex, or 
a geometric mesh, that is, a (pure) 2-dimensional polytopal complex cells. Common cell types 
of a surface mesh are triangles (2-simplices) and quadrangles. Thus a geometric quad mesh is 
a 2-dimensional cubical complex, and a topological one is a cubical 2-dimensional regular CW 
complex. 

In recent years there has been growing interest in volume meshing. Tetrahedral volume meshes 
(simplicial 3-complexes) are well- understood, whereas there are interesting and challenging open 
questions both in theory and practice for hexahedral volume meshes, hexa meshes for short. That 
is, a geometric hexa mesh is a 3-dimensional cubical complex, and a topological hexa mesh is a 
cubical 3-dimensional regular CW complex. 

A challenging open question in this context is whether each cubical quadrilateral geometric 
surface mesh with an even number of quadrangles admits a geometric hexa mesh. In our termi- 
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nology this problem asks whether each cubical PL 2-sphere with an even number of facets admits 
a cubification. Thurston [31] and Mitchell [23] proved independently that every topological quad 
mesh with an even number of quadrangles admits a topological hexa mesh. Furthermore, Epp- 
stein showed in [10] that a linear number of topological cubes suffices, and Bern, Eppstein &; 
Erickson proved the existence of a (pseudo-)shellable topological hexa mesh [5]. 

10.1 Parity change 

Another interesting question deals with the parity of the number of facets of a mesh. For quad 
meshes there are several known parity changing operations, that is, operations that change the 
numbers of facets without changing the boundary. In [5], Bern, Eppstein & Erickson raised the 
following questions: 

(i) Are there geometric quad meshes with geometric hexa meshes of both parities? 

(ii) Is there a parity changing operation for geometric hexa meshes, which would change the 
parity of the number of facets of a cubical 3-ball without changing the boundary? 

Prom the existence of a cubical 4-polytope with odd number of facets we obtain positive answers 
to these questions. 

Corollary 10.1. 

(i) Every combinatorial 3-cube has a cubification with an even number of facets. Furthermore, 
this cubification is regular and even Schlegel. 

(ii) Every combinatorial 3-cube is a facet of a cubical A-polytope with an odd number of facets. 

(iii) There is a parity changing operation for geometric hexa meshes. 

Proof. For (ii) let -F be a combinatorial 3-cube and P a cubical 4-polytope with an odd number 
of facets. By Lemma 3.5 there is a combinatorial 4-cube C that has both F and a projectively 
regular 3-cube G as facets. Let F' be an arbitrary facet of P. Then there is a combinatorial 
4-cube C that has both F' and a projectively regular 3-cube G' as facets. Then the connected 
sum of P and C based on the facet F' yields a cubical 4-polytope P' with an odd number of 
facets, and with a projectively regular 3-cube G" as a facet. The connected sum of P' and G 
glueing the facets G and G" yields a cubical 4-polytope with an odd number of facets, and with 
a projective copy of F as a facet. 

The statements (i) and (iii) follow from (ii) via Schlegel diagrams. □ 

10.2 Flip graph connectivity 

In analogy to the concept of flips for simplicial (pseudo-) manifolds one can define cubical flips 
for quad or hexa meshes; compare [5]. In the meshing terminology the flip graph is defined as 
follows. For any domain with boundary mesh, and a type of mesh to use for that domain, define 
the flip graph to be a graph with (infinitely many) vertices corresponding to possible meshes of 
the domain, and an edge connecting two vertices whenever the corresponding two meshes can 
be transformed into each other by a single flip. 

In this framework, the question concerning a parity changing operation can be phrased as asking 
for a description of the connected components of the flip graph. As an immediate consequence 
of the corollary above we obtain the following result. 

Corollary 10.2. For every geometric hexa mesh the cubical flip graph has at least two connected 
components. 
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11 The next step 



In this paper we are primarily concerned with the reahzation of 2-manifold immersions in terms 
of cubical 4-polytopes, but the higher-dimensional cases are interesting as well. For example, 
one would like to know whether there are cubical 5-polytopes with an odd number of facets. 
(There are no such d-polytopes for d = 6, or for 8 < d < 13; see the [3, Sect. 7].) For this we 
have to realize a normal crossing immersion of 3-manifold into S*^ by a cubical 5-polytope with 
an odd number of quadruple points. Such immersions exist by an abstract result of Freedman 
[13] [1], but more concretely by John Sullivan's observation (personal communication) that there 
are regular sphere eversions of the 2-sphere with exactly one quadruple point [30] [12] and from 
any such one obtains a normal-crossing immersion S-^ with a single quadruple point. 
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